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ABSTRACT 

Vela X-l is the archetype of high-mass X-ray binaries, composed of a neutron star and a massive 
B supergiant. The supergiant is a source of a strong radiatively-driven stellar wind. The neutron 
star sweeps up this wind, and creates a huge amount of X-rays as a result of energy release during 
the process of wind accretion. Here we provide detailed NLTE models of the Vela X-l envelope. We 
study how the X-rays photoionize the wind and destroy the ions responsible for the wind acceleration. 
The resulting decrease of the radiative force explains the observed reduction of the wind terminal 
velocity in a direction to the neutron star. The X-rays create a distinct photoionized region around 
the neutron star filled with a stagnating flow. The existence of such photoionized bubbles is a general 
property of high-mass X-ray binaries. We unveiled a new principle governing these complex objects, 
according to which there is an upper limit to the X-ray luminosity the compact star can have without 
suspending the wind due to inefficient line driving. 

Subject headings: stars: winds, outflows - stars: mass-loss - stars: early-type - hydrodynamics - 
radiative transfer 



1. INTRODUCTION 

A high-mass X-ray binary (HMXB) is a binary star sys- 
tem consisting of a massive luminous hot star (frequently 
OB supergiant) and a compact object, either a neu tron 
star or a black hole (jRemillard feMcClintockl 120061) . A 
fraction of the stellar wind of the luminous hot star is 
trapped in the gravitational well o f the compact object, 
and i s accreted onto its surface ([Davidson fc Ostrikerl 
Il973t lLamers et al.l 119761 ) . Part of the released poten- 
tial energy of accreting material is transformed into X- 
rays, resulting in one of the most powerful stellar X-ray 
sources. Such systems of stars in interaction belong to 
the most valuable astrophysical laboratories. The binary 
nature of the object enables us to determine stellar pa- 
rameters precisely, which subsequently serve as a firm 
base for further study. 

Vela X-l (HD 77581, GPVel) is the archetype 
of high-mass X-ray binaries, consisting of a neu- 
tron s tar and a massive B sup e rgiant dChodil et aTJ 
19671 IBrucato fc Kristianl [1973: IBarziv et al.l 120011: 
Tomsick et al. 1 120101 . The neutron star is a source of 
pulsed X-ray and 7-ray emission wit h a pe riod of 283 s 
(jMcClintock et al.1 119761 : iNorth et al.l 119871 ) , modulated 
both by the orbital motion and stochastic variations 
(|Bildsten et al.lll997ft . The X-rays propagating through 
the hot star wind probe the wind structure, yielding in- 
formation about the m ass-loss rate and the velocity field 
(jWatanabe et al.l 120061 ). The perpetual X -ray variation 
(flaring) reveals t he existence of some structure in th e 
wind - clumping (jDucci et al.ll2009t iFiirst et al.ll2010D . 
On the other hand, X-rays also significantly influence 
the stellar wind, resulting in X-ray photoionization of 
its m aterial (jMacGregor fc Vitellol[l98 2: Wata nabe et al.l 
120061 ) . Because the stellar wind of hot stars is mostly 
driven by the light absorption in the lines of heavier 
elements, the X-ray photoionization may influence the 
wind acceleration. Particularly, the appearance of highly 



charged ions, which absorb the light less effectively than 
low-charged ions, causes the decrease of the radiative 
force. Since this force is responsible for driving the wind, 
the wind flow may subsequently stagnate. 

Numerical studies of stellar winds in HMXBs concen- 
trate mainl y on the multidimens i onal simulation of wind 
accretion felondin et all Il990l: iBlondin fc Wool [19951 : 
iFeldmeier et~aTl 119961 : lHadrava fc Cechural l20ll . while 
the wind driving is simplified using force multipliers that 
take the X-ra y irradiation into an account in an approx - 
imative way (jStevens fc Kallmanl [19901 iStevensI [19911) . 
This is a significant shortcoming, because the X-ray pho- 
toionization affects the radiative force, and consequently 
the amount and velocity of wind material accreted on 
the compact companion. Detailed modelling of ioniza- 
tion and excitation balance in the wind is crucial for the 
understanding of the influence of X-ray photoionization 
on the wind dynamics. The ionization and excitation bal- 
ance should be properly derived using equations of statis- 
tical equilibrium^ To remedy the situation, we provide 
wind models that include the influence of X-ray irradia- 
tion using up-to-date NLTE models. 

2. VELA X-l PRIMARY WIND MODEL 

The applied models of the Vela X-l primary wind are 
based on t he NLTE code with co moving frame (CMF) 
line force (|Krticka fc Kubatl 120101 ). Our models enable 
to selfconsistently predict wind structure just from the 
stellar parameters (the effective temperature, mass, ra- 
dius, and chemical composition). Here we assume that 
the stellar wind of the supergiant component is symmet- 
ric with respect to the binary axis (connecting centers 
of both components) and that the stellar wind in the di- 
rection given by the inclination (f> from the binary axis 

This approach is usually referred to as non-LTE or NLTE and 
it means that the assumption of thermodynamic equilibrium is not 
used for evaluation of the excitation and ionization balance. 
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can be locally described by a spherically symmetric wind 
model. The influence of the neutron star is taken into 
account by its inclusion as the source of external X-ray 
irradiation of the wind. 

Basic parameters of Vela X-l binary system are given 
in Table [T] Binary parameters and the physical pa- 
rameters of binar y members are taken from the spec- 
troscopic analysis (jvan Kerkwiik et~allll995t) . The effec- 
tive te mperature of the supergian t is taken from the ta- 
bles of lStraizys fe Kuriliend (|1981[ ) for the corresponding 
spectral type. The derived value relatively well agrees 
with the determination ba sed on NLTE model atmo- 
spheres dFraser et al.l l2010h . The parameters of the X- 
ray source are a dopte d from observational analysis of 
IWatanabe et al.l ()2006l ) . For ou r calculations we assu me 
the solar chemical composition (jAsplund et al.ll2009( ). 

2.1. Wind model without X-ray irradiation 

For supergiant parameters given in Table [1] we first cal- 
culated NLTE wind model with CMF line force neglect- 
ing the influence of the companion star. The emergent 
surface flux is taken from the H-Hc s pherically sym met- 
ric NLTE model stellar atmospheres (|Kubatll2003l ). The 
predicted wind mass-loss rate 1.5 x 10 -6 M Q year -1 well 
agrees with its estimate derived from the observed X- 
ray spectrum 1.5 — 2 x 10~ 6 Mq year -1 (jWatanab c et al.l 
|2006[) . The predicted terminal velocity 750 kms -1 is 
lower than the observed one 1100 kms -1 (jPrinia et al.l 
Il990f ) . This disagreement may stem either from the 
model simplifications, inaccurate stellar parameters (e.g., 
metallicity), or from the sensitivity of the wind terminal 
velocity to a detailed ioni zation balance in the outer wind 
regions (jPuls et al.ll2000f ). which probably manifests it- 
self in a significant scatter of the ratio of the terminal to 
the escape velocity for the stars with th e same effective 
temperatures (e.g.. iCrowther et al.ll2006D . 

The velocity structure of the wind model without X- 
ray irradiation can be approximated by 
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Note that the representation of the velocity law by 
a polynomial expansion provides better approximation 
than the ordinary velocity law" , b ecause these poly- 
nomi als may form a functional basis (jKrticka fc Kuban 
[201 The exponential term is included for a better fit 
of the velocity law close to the sonic point. The cal- 
culated X-ray opacity per unit mass averaged for radii 
1.5 R* — 5 i?* may be approximated by 
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where A = 10 8 c/^ a x = 0.704 g -1 cm 2 , b 1 = 1.056, 
a 2 = 4.06 x 10 -3 g -1 cm 2 , b 2 = 11.41, and A x = 20.18. 



We stress that A enters as a non-dimensional parameter 
here, which has for a convenience the same value as the 
wavelength in units of A. 

2.2. Modelling of the 2D wind structure 

The full treatment of the problem would essentially 
require a complex s olution of 3D time-dependent hydro- 
dynamic equations (iFriend fc Castorll982l : lBlondin et al.l 
119901 : iFeldmeier et al.lll996[) . However, because the typ- 
ical wind flow-time R^/voa ~ 0.3 day is roughly by a 
factor of 30 shorter than the orbital period, we neglect 
the influence of the orbital motion on the wind structure. 
Consequently, we assume that the stellar wind of the su- 
pergiant component is axisymmetric with respect to the 
binary axis connecting centers of components. Moreover, 
we assume that the stellar wind in the direction given by 
the inclination <p from the binary axis can be modelled by 
a spherically symmetric wind (sec Fig.QJ. Consequently, 
in the orbital plane of the binary our 2D wind model 
consists of sectors of a circle. 

The wind model in each sector is calculated using 
our NLTE wind code. Because the radial wind veloc- 
ity may be non-monotonic in some sectors, we do not 
use the CMF line force for the calculation of these mod- 
els directly, however the line radiative force is given by 
the force cal culated in the Sobolev approximation (e.g., 
lCastorlll974l) multiplied by the ratio of the CMF to the 
Sobolev line force derived from the wind model that ne- 
glects the radiation from the companion star (Sect. I2~l"j) . 

The influence of the neutron star is taken into account 
only by inclusion of its X-ray radiation due to the wind 
accretion on the neutron star. This radiation irradiates 
the supergiant and interacts with its wind. To describe 
this effect, we add an additional term J x (v,r,d) to the 
specific intensity J(v, r, d) in the form 



J x (v r d) = I^L C -T(v,r,d) 
J 1 ' ' > ~ 16^2 6 



(4) 



where the optical depth along the given ray is (z mea- 
sures the distance along this ray) 



T(v,r,d) = / k(v, z)p(z)dz, 
Jo 



(5) 



L x {y) is the X-ray luminosity per unit of frequency, d is 
the distance from the given point in the supergiant wind 
region to the surface of the neutron star, z) is the ab- 
sorption coefficient per unit of mass, and p(z) is the wind 
density. The distribution of em ergent X-rays L x (v) i s 
approximated by the power law (|Watanabe et al.ir2006| ). 
Energies higher than 3 keV, which are well above the ion- 
ization energies of all included ions were not considered 
in the model. The absorption coefficient and the den- 
sity in Eq. ([5]) can be derived from models for individual 
sectors. However, to simplify our approach, for the cal- 
culation of J x {v) we use fits following from Eqs. ([1]), © 
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Table 1 

Parameters of Vela X-l (HD 77581) binary system 



Parameter 


Value 


Reference 


Binary 


Separation D 


53.4 R Q 


van Kcrkwiik ct al. (1995) 


Period P 


8.96 day 


van Kcrkwiik ct al. (1995) 


Supergiant 


Spectral type 


B0.5Iae 




Radius R* 


30 Rq 


van Kerkwiik et al. (1995) 


Mass M 


23.5 M 


van Kerkwiik et al. (1995) 


Effective temperature T e g 


27 000 K 


Straizvs &; Kurilicnc (1981) 


Wind mass-loss rate M 


1.5 X 10" 6 M Q year- 1 


this work 


Wind terminal velocity Voo 


750 km s" 1 


this work 


Neutron star 


Mass 


1.88 M 


van Kcrkwiik ct al. (1995) 



X-ray luminosity L x 3.5 X 10 36 ergs -1 Wa.tana.be ct al. (2006) 





-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

xl D 

Figure 1. The optical depth for 830 cV between the neutron star 
and a given point plotted in the plane containing the binary axis. 
The distances are given in the units of binary separation D with a 
small ruler plotted in units of supergiant radius R* . Note that the 
optical depth plotted in the figure is only an approximative one, 
because it was calculated using fits Eqs. JTJ, J3]l to wind solution 
without any external X-ray irradiation. Anyway, the main trends 
are the same as in a more precise 2D model. 

where the relation between the distance along the ray 
z and the radius r is derived from the geometry of the 
problem, and v{4>) is an average velocity in the velocity 
plateau which occurs due to X-rays (see Fig. [3]). 

The models describing the wind in different sectors 
with inclination (f> with respect to the binary axis are 
calculated for a sequence in <fi with a step of 10°. 

3. WIND STRUCTURE IN 2D 

The X-ray source located on the surface of the neu- 
tron star influences the ionization state of the supergiant 
wind. The influence is stronger if the X-ray optical depth 
between a given point in the wind and the neutron star 
surface is lower. To illustrate this, in Fig. Q] we plot the 
optical depth (Eq. ([5])) in a plane containing the binary 
axis. Due to the assumed symmetry of the problem, the 
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Figure 2. Influence of X-ray photoionization on the ionization 
state of selected elements. Calculated for a sector with = 50° . 

optical depth is axisymmetric with respect to the binary 
axis. The X-rays strongly penetrate the wind that di- 
rectly faces the neutron star. On the other hand, due 
to geometrical reasons, the radial wind streams that are 
significantly inclined with respect to the neutron star are 
affected by X-ray radiation at larger distances from the 
primary. 

As a result of the X-ray photoionization of the wind, 
lower ionic states are effectively destroyed and higher 
ionic states appear in a nonncgligiblc amount. This can 
be seen from Fig. [3J where we compare the ionization 
fraction of selected ions in the model with and without 
X-ray irradiation. The X-rays influence the wind ioniza- 
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Figure 3. The radial wind velocity for different inclinations <j> 
with respect to a neutron star. The position of the neutron star is 
marked in the graph. 

tion state in the region where the X-ray optical depth 
between a given point and the neutron star is low, r < 1, 
i.e, to a distance comparable to an orbital separation D. 
X-rays are not able to penetrate the wind close to the 
supergiant star (since r 3> 1 if we aim to approach the 
supergiant surface), consequently the ionization state of 
material there is not affected by X-rays. 

Because the stellar wind of hot stars is accelerated 
by the light absorption in the lines of heavier element 
ions, any change of the wind ionization state affects 
the accelerating radiative force. On average, ions with 
higher charge are less effective in driving the wind, ba- 
sically due to the lower number of their spectral lines. 
For a weak X-ray irradiation the emergence of new 
ionization stat es causes a slight incre ase of the radia- 
tive force (cf. iKrticka fc Kubatl 120090 . On the other 
hand, for a strong irradiation when the degree of ion- 
ization is higher the radiative for ce significantly de- 
creases (jMacGregor fc Vitellol QJ982) . Because this ef- 
fect appears in the outer wind region and the wind 
mass-loss rate is determined in the inner wind region, 
only the wind velocity (see Fig. |3l). and not the wind 
mass-loss rate (jMacGregor fc Vitelld[l982l ) is affected by 
X-ray photoionization. This is especially apparent for 
lower wind inclinations (<f> < 40°), where the wind ra- 
dial velocity becomes nonmonotonic. The decrease in 
the radiative force is so strong there that the acceler- 
ating wind solution is not possible any more, and the 
wind velocity switches to decelera ting overloaded solu- 
tions (jFeldmeier fc Shlosmanl 120001 ) with a typical kink 
in the velocity profile. This has further consequences for 
even lower wind inclinations. The radiative force in the 
rays that closely pass the neutron star {<t> < f5°) is so 
strongly affected by the X-ray photoionization that the 
radiatively driven wind in these directions never reaches 



the neutron star. 

As the result of X-ray photoionization, a region (bub- 
ble) around the neutron star is created, in which the 
wind stagnates at low velocities, and consequently high 
densities. These results explain the detection of low- 
sp eed wind regions in the direction of the neutron star 
bv lWatanabe et al.l (|2006f ). 

4. DISCUSSION 
4.1. X-ray luminosity 

Part of the mass which leaves the primary star via its 
wind can be accreted onto the compact component. Its 
potential energy is then released in the form of X-ray ra- 
diation. Let us adopt an approximative Hoyle-Lyttleton 
treatment to estimate the X-ray luminosity. For moving 
medium, the radius from which the matter can be ac- 
creted onto a point mass is given by the relation thl = 
2GMx/v 2 , where Mx is the mass of the compact compo- 
nent (neutron star), v 2 = v^ind + v^ lh , Vorh * s ^ e or bital 
velocity of the compact component, and w W i n d is the wind 
velocity at the distance of the compact component (de- 
noted by D). Given u or b = 300 km s -1 and our simulated 
wind structure, we find that the neutron star is able to 
collect matter only from a narrow cone defined by the 
value of 4> < 15 (for <\> = 15 and the distance D we have 
«wind = 32.0 kms -1 and t^hl = 0.17D). We can express 
the relation between the accretion rate M acc and the 
mass loss rate M from the supergiant as (jWatanabe et al.1 
120061 ) Af acc = Afr| L /(4£> 2 ). Then the X-ray luminosity 

L x = GM x M acc /R x is L x = 8.8 x lO^ergs" 1 . This 
value is approximately one order of magnit ude higher 
than the observed one ijWatanabe et all [20061 ) . The rea- 
son of this difference is the fact that only a small fraction 
of matter in the accretion cone gets finally accreted and 
the rest falls back to the surface of the primary as a con- 
sequence of the flow stagnation between the supergiant 
and the neutron star. 

4.2. Existence of two types of solutions 

The wind equations allow the existence of two types 
of solutions giving different X-ray luminosities and wind 
velocities. The solution presented here appears in the 
case of strong X-ray source, which significantly affects the 
wind ionization state and consequently also the radiative 
force. This results in a slow wind that can be accreted 
by the neutron star in large amounts producing a strong 
X-ray source. 

Another type of solution may occur in the case of a 
weak X-ray source that does not significantly influence 
the wind ionization state. The radiation force becomes 
higher, similar to that without any X-ray irradiation. 
This results in faster outflow roughly corresponding to 
the "no X-rays" case in Fig. [3] Due to the dependence 
of the accretion rate on the velocity via M acc ~ v -4 and 
having roughly two times higher v, the X-ray luminosity 
is by an order of magnitude lower in this case (since Lx ~ 
Mg.cc). This agrees with adopted assumption of weak X- 
ray source. Consequently, also this is a possible solution 
of wind equations. 

Different types of solutions may appear in different sys- 
tems, or even an external perturbation may cause switch- 
ing between these two wind solutions in a particular bi- 
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Figure 4. The location of forbidden X-ray lumi nosit ies as a func- 
tion of the optical depth parameter (Eqs. (5), fllPI) ), Observed 
parameters of individual HMXBs (crosses) are located outside the 
highlighted forbidden area. Here w e plotted points (denot ed by 
number s) for binaries V ela X-l (1, [Watanabe ct al. 2006j ) , Cy g 

- ' Prat eTaLI 2003) , 
Ribo et al. 2006 . 



X-l (2. lYan et all 12001). ICR J19140+0951 (3, 
Wray 977 (4. IKaper et al.ll200l), 4U 2206+54 (5,L 
4U 0114+6 5 (6. IReig et al.l 119961; IMukhcricc & Paul 200 ^ 4U 
1907+09 (7. ICox c t al. 2003), and LS 5039 (8, ICasares et al.|[2o05T) . 
The plot of the maximum luminosity was derived using the opti- 
cal depth between the neutron star and the wind critical point 
rx (Vela X-l) = 9. 

nary system, possibly contributing to the variability of 
X-ray luminosity. This variability might be accompa- 
nied by the variability of the distribution of emitted X- 
rays, if the accretion regime changes with accretion rate 
(jShakura et alj|2012j) . 

4.3. Implication for other HMXBs 

For a slightly lower mass-loss rate or for a slightly 
higher X-ray luminosity than assumed here the X-rays 
could penetrate deeply into the stellar wind and signifi- 
cantly influence the ionization state at the wind base re- 
sulting in the decrease of the radiative force. This could 
lead to the disruption of the stellar wind and a significant 
decrease of X-ray luminosity, possibly providing another 
contributi on to the overall X-ray va riability observed in 
Vela X-l (jKrevkenbohm et al.ll2008l) . 

This means that there exists a maximum X-ray lumi- 
nosity the compact star can have for a given geometry of 
the system and wind mass-loss rate. Assuming that the 
Vela X-l luminosity is the maximum one, the maximum 
allowed X-ray luminosities for other HMXBs L™ ax can be 
derived using the optical depth (compare with the ion- 
ization parameter f. lTarter et al.lll969l ) rx = J R npdr ~ 

M(D — i?*)/(w 00 i5i?*) and corresponding parameters of 
Vela X-l as 



- x ( - = L x (Vela X-l)e 
or, in scaled quantities 



-T X (Vela X-l) 



M 



(9) 



log(L x lax /l erg) = 32.6+3.9 

i / (DR*)/(D - R 



WsOkms" 1 / 



1.5 x lO^Moyear- 1 I 

• (io) 



68.5 R ( 



Here M, Voa, i?*, and D are the wind mass- loss rate, the 
terminal velocity, radius of luminous component, and the 
binary separation for individual HMXBs. 



The results of observations for other HMXBs collected 
in Fig. |4] clearly support the picture that there exists a 
maximum allowed X-ray luminosity, which depends on 
the wind and geometry parameters. Some systems lie 
close to the boundary of the forbidden area, whereas 
others have typically higher wind mass-loss rate M > 
5 x lO _6 M0year _1 that keeps them out of the bound- 
ary. The position of individual stars in this diagram may 
vary with time due to the existence of two possible solu- 
tions of wind equations (as discussed in Sect. I4.2|) . 

4.4. Mass-loss rate determination 

Our results support the mass-loss rate predic- 
tions based on up-to-date wind models in two ways. 
First, our mass-loss rate prediction of HD 77581, 
1.5 x 10 -6 Mq year -1 , agrees with the value estimated 
from X-ray spectrosco py, 1.5 — 2 x 10~ 6 MQyear -1 
([Watanabe et al.l 120061 ) . Moreover the wind mass- loss 
rate cannot be lower than this value, because a decrease 
of the wind mass-loss rate would lead to lower X-ray 
opacity, stronger wind X-ray photoionization close to the 
star, and even more significant reduction of the radiative 
force. This would finally cause a disruption of the wind 
and a disappearance of X-ray emission. This imposes a 
strong lower limit on the observational wind mass-loss 
rate estimates, that is in agreement with current mass- 
loss r ate predictions (|Vink et alj|2001t iKrticka fc Kubatl 
I20"l0l) . 

4.5. Limitations of the present models 

The consistent inclusion of the X-ray irradiation, which 
is an advantage of our models, also determines their 
shortcomings. The coupled solution of NLTE and ra- 
diative transfer equations is significantly time-consuming 
even in ID. Its inclusion into multidimensional time- 
dependent simulations is likely beyond the possibilities 
of current computers. Consequently, while one part of 
the problem solution is treated in detail, the second one 
is simplified. The correct picture of the flow in the 
HMXBs should take into an account the results of both 
approaches: the stagnation of the flow in the direction 
towards the neutron star, which is studied in this paper, 
and a comple x 2D picture of the wind accretion on the 
neutron star (IBlondin et alJll990t iBlondin fc Wodll995l: 
iFeldmeier et al.lll996| ).~ 

On the other hand, there are effects that are not de- 
scribed by any of the available models. It is well es- 
tablished that the hot star wind is inhomogeneo us on 
small scales (clumped, see lHamann et al.l [2008). In 
the case of HMXBs, clumping (which favors recombi- 
nation) may affect th e region in which th e photoion- 
ized bubble is formed (Oskinov a et al.ll2012l). However , 
the hydrodynamical simu lations (jFeldmeier 
iRunacres k, Owockl l2002f) predict that clumping starts 
above the critical point of the wind solution, conse- 
quently, it does not affect wind mass-loss rates and ter- 
minal velocities. Thus we expect that clumping does not 
significantly influence the results of our models. 

The wind inhomogeneities likel y cause the high 
variability of the X-ray source (|Furst et al.l 120101 : 
lOskinova et al.ll2012D . The calculation of the wind ion- 
ization should in fact account for such time-dependent X- 
ray photoionization. However, because the typical flow 
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time of the wind is longer than the typical timescale of 
X-ray variability, we expect that our models are able 
to reproduce the mean effect of the X-ray photoioniza- 
tion. On the other hand, the variable X-ray ionization 
source causes temporal changes of the radiative acceler- 
ation, providing external perturbation which may con- 
tribute to the natural wind clumping. As a result, the 
primary wind may be more clumpy than the wind of a 
similar single supergiant. 

5. CONCLUSIONS 

We provide detailed numerical models of the influence 
of X-rays on the supergiant wind in the Vela X-l binary 
system. The effect of X-ray photoionization on the ra- 
diative force and wind dynamics has never been studied 
using appropriate NLTE wind models. The X-rays pho- 
toionizc the wind and destroy the ions responsible for 
the wind acceleration. This results in flow stagnation 
in the vicinity of the neutron star, which was identified 
in observations. For a sufficiently strong X-ray source 
the wind that directly faces the neutron star falls back 
on the mass losing star and never reaches the compact 
companion. 

We have shown that there is an upper limit to the 
X-ray luminosity the compact star can have without dis- 
rupting the stellar wind. For a higher luminosity than 
the limiting one the decrease of the wind acceleration is 
so strong that no wind material would reach the neutron 
star. This theoretical picture of the maximum X-ray lu- 
minosity is supported by observation of many high-mass 
X-ray binaries. 

The wind equations allow the existence of two types 
of solutions with different X-ray luminosities and wind 
velocities. The case of a strong X-ray source, which sig- 
nificantly affects the wind ionization, leads to accretion 
of slow wind in large amounts resulting in a strong X- 
ray source. On the other hand, a weak X-ray source that 
does not significantly influence the wind ionization re- 
sults in accretion of fast wind in low amounts producing 
a weak X-ray source. Different types of solution may ap- 
pear in different binary systems, or perturbations may 
cause switching between these two types of wind solu- 
tions contributing to the X-ray variability. 

The predicted mass-loss rate agrees with the value es- 
timated from X-ray spectroscopy. Moreover, the wind 
mass-loss rate cannot be lower than this value, because 
a decrease of the wind mass-loss rate would lead to the 
disruption of the wind and disappearance of the X-ray 
emission. This supports the reliability of current mass- 
loss rate predictions. 
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